Current drag in capacitevly coupled Luttinger constrictions 
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We study the current drag in the system of two electrostatically coupled finite ID electron 
channels. We present the perturbation theory results along with the results for two non-perturbative 
regimes. It is shown that the drag may become absolute, that is, the currents in the channels are 
equal in a finite window of the bias voltages. 



The recent experiments with two capacitevly coupled 2DEG demonstrated that, owing to Coulomb interaction, 
electrons moving in one of the 2DEG drag the electrons of the adjacent 2DEG. Therefore the system work as a d.c. 
current transformer. The transformation coefficient was however fairly small. Later it has been shown that in a 
specially designed coupled ID arrays of ultrasmall tunnel junctions Coulomb interaction may lead to the absolute 
current drag This means that the electric currents in two capacitevily coupled circuits are equal in magnitude in 
a certain region of voltages applied to the circuits, and the system may work as a current copier. This prediction has 
been very recently confirmed experimentally. 

The mechanism of the absolute current drag in small tunnel junctions is however quite different from the momentum 
transfer physics that is responsible for the current drag in the 2DEG systems. The aim of the present work is to show 
that in the ID electron channels, the momentum-transfer mechanism can also lead to the absolute current drag. We 
consider a system of two channels of finite length L that are adiabatically connected to the reservoirs of effectively non- 
interacting electrons (Fig.l). The channels are assumed to be coupled only by the Coulomb interaction. Further, we 
assume no impurities in the channels that would cause electron backscattering. Being inspired by recent advances Q| 
in fabrication of relatively long one-channel quantum wires, we hope that such systems can be successfully fabricated 
in the nearest future. 

To get an intuitive feeling of how the absolute current drag may occur in such a system, let us consider the case of 
strong electron repulsion. In this case, electrons in each channel form a rigid Wigner lattice. Provided the channels 
are close to each other, the repulsion coordinates not only the positions of electrons inside the same channel but also 
the positions of neighboring electrons in the other channel. (Fig. 1) Now, if the electrons in one of the channels move, 
electrons in the other channel must follow their motion, so that the electric currents in the channels are equal. (Fig. 



Such a simple model says very little about realistic ID channels where electrons are subjected to strong quantum 
fluctuations. Below we develop a consistent theory of current drag which accounts for these fluctuations. The main 
prediction of the theory is that the (almost) absolute current drag survives quantum fluctuations. It can occur even 
if the repulsive interaction is weak provided the channels are sufficiently long. 

We model each of the two coupled ID conductors as "Luttinger constriction" with spinlcss electrons, i.e., use the 
standard Hamiltonian 



Here Q{x) is the displacement of electrons in the conductor normalized in such a way that the local fluctuation of 
electron density is 5n = — VG/tt, and n(a;) is the conjugated momentum density: [6(x), n(a:')] = mh6{x — x'). The 
velocity u of excitations in the conductors depends on x: in the constriction region of length L the electron-electron 
interaction gives an essential contribution to compressibility and u is larger than outside of the constriction, where 
interaction is screened out and u is equal to the Fermi velocity up- Parameter g in the Hamiltonian (|^) characterizes 
the strength of the interaction, g = uf/u. 

The relevant part of the interaction of the two conductors responsible for the momentum transfer between them 
can be written as 
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where subscripts 1,2 refer to the two conductors, kpj = nrij are the Fermi wavevector, and J is the interaction 
constant, precise value of which depends on the specific electrostatic configuration of the system. In particular, it 
depends on the distance d between the two channels, decreasing roughly as exp(— dn) at large d, where n is the average 
electron concentration in the channels. 

When the coupling is weak, electrons in the two channel move independently and can carry generically different 
currents Ii and l2- When the currents are different however, momentum transfer between the channels leads to the 
appearance of an extra voltage drop 6Vj across each channel, in addition to the "ideal" Landauer value hlj/e^ . The 
extra voltages 5Vj can be found from the difference of chemical potentials at the ends of the channels: 

5Vj = ^(5nj(-L/2) - 5jIj{L/2)) . (3) 

where k's are compressibilities of electron liquids in the channels. Writing the equations of motion for phases that 
follow from the Hamiltonian H = Hi + H2 + H^- 

Qj - [u]Q'j\' = {-iyjuFsm[2[kFi - kF2)x + 2[Qi{x) ~ Q2{x))] , (4) 
and using the fact that in the stationary case (G)) = 0, we get that 5V2 — —5Vi = V , and: 

V=(' dx^^, F{x)^ J{sin{2{kFi~kF2)x + 2{Qi{x)~Q2{x))) . (5) 
J-L/2 e 

In the first order of the perturbation theory in He we get: 
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where k = 2(fci?i — kF2), and A{x,x',t) — ([1^, (^']+ — ip"^ — ip''^)/2, B{x,x',t) — {[(p,ip'])/2 . Here (p = ip{x,t), 
if' = ip{x',0), and the phase ip{x,t) is introduced by the relation: 

2ieiix,t)-e2{x,t)) = flt + if{x,t), n = 2^{h~l2)/e. (7) 

The average (...) in the definitions of the correlators A and B is taken over the equilibrium fluctuations of ip. The 
assumption of equilibrium is legitimate since after we separated the linearly growing part VLt in eq. (^) the dynamics 
of 0's in the zero order in He does not have any perturbations that would drive the system out of the equilibrium. 
The equilibrium correlators A, B can be expressed in a standard way in terms of the spectral densities Pj{x, x' ,uj): 

A{x,x',t)±B{x,x',t)^2^ duj e'-'^'^Pj{x,x',uj){coih-^±l) - Pj{x,x,oj) - Pj{x',x',uj) , (8) 

where spectral densities are pj — (l/7r)ImG^'' (x, a;', w). The retarded Green's functions G^"* of the phases &j satisfy 
the equations (^) with J = 0, radiative boundary conditions at a; — > ±00, and delta-functional source term in the 
right-hand-side. Such equations can be solved explicitly for our model of piecewise constant u{x) and give for x, x' 
inside the constriction: 

, , (1 + (7^) cos;^(a; — a;') + (1 — 0^) cosi^i cosi'fa; -I- a:') 

p(x,x'uj)^—- — ^ — ^ -, v = uju. (9) 

' ' ^ uj (1 + 52)2 sin2i.i + 452 cos2i.L ' ' ' 

Equation (^) combined with eqs. (|^) and allows us to find the induced voltage ^ as a function of the current 
difference / = /i — /2 in various regimes. We consider first identical channels at zero temperature. In the non- 
interacting case (7 = 1 we get then from eq. (^): 

■kPLuVL^ f^^/"^ dx smx - xcosx 



^-^^/. x^ ^'-m^^ ' (^0) 



where Ef is the cutoff frequency of the excitation spectrum of the conductors, being of the order of Fermi energy. 
Equation (|l^) shows that at small current differences il, for which Vl <C u/L, the voltage V grows as Vt^, and this 
growth slows down to il^ at large Vl. 
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At small frequencies the spectral density is dominated by the leads and approaches its non-interacting value 
l/2w at arbitrary interaction strength inside the constriction. Since at small current differences fl the system does 
not have enough energy to excite high-frequency modes in the constrictions this means that, the dependence V (x ft^ 
is valid at sufficiently small f2 for any g. However the range of the validity of this dependence shrinks with growing 
interaction strength (decreasing 17) to f2 <C gu/L. 

To see how V behaves at large current differences {il ^ u/L) we can average the rapid oscillations of spectral 
density with the period u/L and also neglect the terms that oscillate rapidly with x,x'. The spectral densities 
reduce then to pj{x, x' ,uj) = gj cosv^x — x')/2uj, and we get for V: 

ig-2 
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We see that the interaction inside the constriction changes drastically the behavior of y as a function of giving 
rise to power-like dependence. The voltage V decreases at large I ioi g < 1/2, so that there is optimum / ~ eu/L 
for which the induced voltage is maximum. If we take oscillating terms into account, the decrease at large current 
differences is non-monotonic, the voltage oscillates with f2 with the period u/L. This resembles density of states 
behavior in finite Luttinger constriction ||]. All these features can be seen in Fig. 2 which presents the results of 
numerical calculations of 1^ as a function of Ii — I2 from Eqs. 6-9. 

Another important feature of Eq. |l| is that the voltage V grows with increasing L at any given current difference 
fl regardless of g. This suggests that Eq. ^ is valid only for sufficiently short channels, when the voltage V is small 
perturbation. In longer channels the distribution of the voltage along the channel and associated with it difference of 
electron densities in the two channels should be determined self-consistently. The two quantities are related through 
the differential version of Eq. ^, which for identical channels has the form V'{x) = —Kk' {x)/2ne. The variation of the 
density difference can be described by the steady-state form of Eq. ^ which can be written as follows: 

k'ix) = ^F{k{x),n). (12) 
u 

Here the force F is given by Eq. |ll| but now with local k{x). In simple terms, Eq. |l^ describes how the friction 
between the moving liquids compresses one of them and stretches another one, and how this compression/stretch 
affects the magnitude of the friction itself. 

Equation with the force given by eq. pi] shows that in long channels, L ^ Ln, where L'^^ = l&g,Pu{Vl/2EF)'^^ 
is characteristic length of variation of k{x), the density difference approaches the value fl/iru and saturates. This 
means that the extra voltage difference drops at distances of the order of Lq near the edge of the interacting region 
and is equal to h{Ii — l2)/e^2g. The total voltage drops across the channels are given then by "quasi Landauer" 
relations 

e^Vi/h^h + ^ih-h), e^V2/h = l2 + ^{l2-h)- (13) 
2g Zg 

We see from this equations that the transconductance of the two channels depends strongly on the interaction strength 
g inside the channels, and for strong interaction, g 0, the currents in the two channels are almost equal in the wide 
range of bias voltages. 

Another nonperturbative regime in which the currents in the two channels are almost precisely equal is the limit 
of strong coupling of the two short channels with L ^ Lq . This regime corresponds to the classical limit of the sine- 
Gordon model. In the limit of g ^ 0, one can treat cp as (almost) classical variable that takes on 27rxinteger values. 
In our original model, this means that the Luttinger liquids, or electron chains, are strictly correlated. However, 
they can move together and carry equal currents in both constrictions. The relation between current and voltages in 
this case is readily obtained from the fact that, as follows from Eq. (||), the deviations from Landauer voltages are 
opposite in both channels, SVi — —6V2. From this we obtain Ii = I2 — e^iVi -f V2)/2h. 

This is, of course, an ideal situation, and at finite g and L the currents are not precisely equal. Let us show 
that, in the wide region of parameters, the current difference is exponentially small. To this end, we switch from the 
Hamiltonian to Lagrangian in imaginary time. This representation is standard for calculations of rates of classically 
forbidden processes. The Lagrangian describing the relative motion of the two electron liquids within the interaction 
region takes the standard sine-Gordon form. 



S[Lp{x,t)] = 



dxdt 
2-K 
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Equation [7| shows that the shift of 4> at the boundary by 47r corresponds to transfer of one electron into the channel 1 
accompanied by the transfer of an electron out of the channel 2. This suggests that, in order to describe the interaction 
of the electron channels with the electrodes, we should add to the Lagrangian the boundary term 

= j ^ [(^^^ - ^'l)vi-L|2) + {^,[ - , (15) 

where fii'^ are the chemical potentials of the left/right electrode of the channel 1 and 2. 

What are the processes that spoil the ideal current copying? Obviously, those are slips of (p by 27r. If the voltages 
applied to the constriction are sufficiently low, the slip should happen at once throughout the interacting region. 
Otherwise a soliton would have been created within the interacting region as a result of the slip. This would cost 
energy Es — Ju/irg and is forbidden by the energy considerations. The sudden slip throughout the interacting 
region may be viewed as a virtual soliton transfer via the " barrier" Es provided the barrier length exceeds the soliton 
size ^Jujjg. The rate of such a process is determined by the saddle-point configuration of 0(a;, r) as sketched in Fig. 
3. The action is contributed by two soliton "walls" and the boundary term (/i = max((^f' — )/2, (/if — /i2^)/2)), 

S = IEs^Jt"^ + {Ljuf - 2Tfi. (16) 

Minimization of this expression with respect to r yields 

So^2^E^~ti^L/u, (17) 

and h — h oc cxp(— 6*0). We see that the current difference is suppressed exponentially provided L » u/Eg and 
Mi.-R < ^s- The common wisdom of the sine-Gordon model suggests that these results are in fact valid not only in 
the classical limit g ^ Q but also at any 5 < 1, provided Es is properly renormalized corresponding to the actual 
soliton energy Es{g). The concrete formula for Es{g) reads 

Es{g)^Es{EslEpy''~^ (18) 

with the coefficient depending upon cutoff details. There is an interesting analogy between the soliton and Dirac 
particle with the dispersion relation E{k) — ^JE^+v?}? in the interacting region and E{k) = uk beyond that. At 
energies E < Eg the particle would tunnel through the interacting region. The energy- and L-dependence of such a 
tunneling process coincide exactly with that of Eq. |l^. 

In conclusion, we have considered current drag in the system of two capacitevely coupled Luttinger constrictions. 
We have presented perturbative results and have found two nonperturbative regimes. In one of the regimes, the 
system works as an ideal current copier of exponential accuracy. 

We wish to thank A. A. Odintsov for his contribution at preliminary stages of this work. We are indebted to S. 
Tarucha for his communications that gave an impetus to the present research, G.E.W. Bauer, M.P.A. Fisher for very 
instructive discussions of the results. This work has been made possible by financial support of North Atlantic Treaty 
Organization, grant No. 950279. 



[1] P.M. Solomon, P.J. Price, D.J. Frank, and D.C. La Tulipe, Phys. Rev. Lett. 63, 2508 (1989); T.J. Gramila, J.P Eisenstein, 

A.H. MacDonald, L.N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 66, 1216 (1991). 
[2] D.V. Averin, A.N. Korotkov, and Yu.V. Nazarov, Phys. Rev. Lett. 66, 2818 (1991). 

[3] M. Matters, J.J. Versluys, and J.E. Mooij, Phys. Rev. Lett. 78, 2469 (1997);P. Delsing, D. B. Haviland, and P. Davidsson, 

Chechoslovak journal of physics, 46, 2359 (1996). 
[4] S. Tarucha, T. Honda, and T. Saku, Sol. St. Comm. 94, 413 (1995); A. Yacoby, H.L. Stormer, and K.W. West, Phys. Rev. 

Lett. 77, 4612 (1996). 

[5] F.D.M. Haldane, Phys. Rev. Lett. 47, 1840 (1981), see for review J. Voit, Rep. Prog. Phys. 58, 997 (1995). 
[6] Yu.V. Nazarov, A.A. Odintsov, and D.V. Averin, Europhys. Lett. 37, 213 (1997). 
[7] S. Coleman, Phys. Rev. D 11, 2088 (1975), and references therein. 

[8] R.F. Dashen, B. Hasslacher, and A. Neveu, Phys. Rev. D 11, 3424 (1975); A. Luther, Phys. Rev. B 14, 2153 (1975). 

FIG. 1. The system under consideration. Two ID channels are open to reservoirs and separated by inslulating barrier. 
The electrostatic repulsion coordinates the positions of electrons (black circles) in the channels, that may lead to the absolute 
current drag. 
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FIG. 2. The perturbative regime. The induced voltage difference is normalized by Vo{g) = J^uLEp^{u/EfL)*^ ^.The 
curves a,b,c,d,e correspond to the values oi g = 0.1, 0.3, 0.5, 0.6, 1.0 respectively. 



FIG. 3. Soliton tunneling via the interacting region. The optimal configuration of i^(a;,r) correspond to three regions of 
almost constant </? separated by V-shaped soliton "wall" (thick line). 



5 



